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Non-perturbative Hamiltonian light-front quantum field theory presents opportunities and challenges that 
bridge particle physics and nuclear physics. Fundamental theories, such as Quantum Chromodynmamics (QCD) 
and Quantum Electrodynamics (QED) offer the promise of great predictive power spanning phenomena on all 
scales from the microscopic to cosmic scales, but new tools that do not rely exclusively on perturbation theory 
are required to make connection from one scale to the next. We outline recent theoretical and computational 
progress to build these bridges and provide illustrative results for nuclear structure and quantum field theory. As 
our framework we choose light-front gauge and a basis function representation with two-dimensional harmonic os- 
cillator basis for transverse modes that corresponds with eigensolutions of the soft- wall AdS / QCD model obtained 
from light-front holography. 



1. Introduction 

Major goals of Hamiltonian light-front field 
theory include predicting both the masses and 
transitions rates of the hadrons and their struc- 
tures as seen in high-momentum transfer exper- 
iments. So far, the approaches in a discretized 
momentum basis [ [1] and in transverse lattice [ 
[21 Oil] have shown significant promise. The main 
advantages of Hamiltonian light-front quantum 
field theory are: (1) one evaluates experimental 
observables that are non-perturbative and rela- 
tivistically invariant such as masses, form factors, 
structure functions, etc.; (2) one evaluates these 
quantities in Minkowski space; and (3) there is no 
fermion doubling problem. 
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Here we present a basis-function approach 
that exploits recent advances in solving the non- 
relativistic strongly interacting nuclear many- 
body problem [[SHU]. Both light-front field theory 
and nuclear many-body theory face common is- 
sues within the Hamiltonian approach - i.e. how 
to (1) define the Hamiltonian; (2) renormalize to 
a finite space; (3) solve for non-perturbative ob- 
servables while preserving as many symmetries 
as possible; and, (4) take the continuum limit. 
Nevertheless, Ken Wilson has assessed the ad- 
vantages of adopting advances in quantum many- 
body theory and has long advocated adoption of 
basis function methods as an alternative to the 
lattice gauge approach [\7\. 
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2. Ab initio Hamiltonian approaches to 
quantum many-body systems 

To solve for the properties of nuclei, self-bound 
strongly interacting systems, with realistic Hamil- 
tonians, one faces immense analytical and com- 
putational challenges. Recently, ab initio meth- 
ods have been developed that preserve all the 
underlying symmetries and converge to the ex- 
act result. The basis function approaches that 
we adopt here, No Core Shell Model (NCSM) [[5] 
and No Core Full Configuration (NCFC) meth- 
ods [H>], are among the several methods shown to 
be successful. The former adopts a finite basis- 
space renormalization method and applies it to 
realistic nucleon-nucleon (NN) and three-nucleon 
(NNN) interactions (derived from chiral effective 
field theory) to solve nuclei with Atomic Num- 
bers A = 10—13 [ [8]. Experimental binding 
energies, spectra, electromagnetic moments and 
transition rates are well-reproduced. The lat- 
ter adopts a realistic NN interaction that is suf- 
ficiently soft that renormalization is not neces- 
sary and binding energies obtained from a se- 
quence of finite matrix solutions may be extrap- 
olated to the infinite matrix limit. Owing to uni- 
form convergence and the variational principle, 
one is also able to assess the theoretical uncer- 
tainties in the extrapolated result. One again ob- 
tains good agreement with experiment. The pri- 
mary advantages of these methods are the flexi- 
bility for choosing the Hamiltonian, the method 
of renormalization/regularization and the basis 
space. These advantages direct us to adopt the 
basis function approach in light-front quantum 
field theory. In a NCSM or NCFC application, 
one adopts a 3-D harmonic oscillator for all the 
particles in the nucleus (with harmonic oscillator 
energy fi), treats the neutrons and protons in- 
dependently, and generates a many-fermion basis 
space that includes the lowest oscillator configu- 
rations as well as all those generated by allowing 
up to N max oscillator quanta of excitations. The 
single-particle states are formed by coupling the 
orbital angular momentum to the spin forming 
the total angular momentum j and total angu- 
lar momentum projection Mj. The many-fermion 
basis consists of states where particles occupy the 
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Figure 1. Calculated ground state energy of 12 C 
for N max — 2—8 (discrete points) at selected val- 
ues of the oscillator energy, CI. For each CI, the re- 
sults are fit to an exponential plus a constant, the 
asymptote, which is constrained to be the same 
for each curve [|6j. We display the experimental 
ground state energy and the common asymptote. 
Figure from Ref . [ [TO] . 



allowed orbits subject to the additional constraint 
that the total angular momentum projection Mj 
is a pre-selected value. For the NCSM one also 
selects a renormalization scheme linked to the 
many-body basis space truncation while in the 
NCFC the renormalization is either absent or of 
a type that retains the infinite matrix problem. 
Non-perturbative renormalization has been de- 
veloped to accompany these basis-space methods 
that preserve all the symmetries of the underlying 
Hamiltonian including highly precise treatments 
of the center-of-mass motion. Several schemes 
have emerged with impressive successes and cur- 
rent research focuses on detailed understanding of 
the scheme-dependence of convergence rates (dif- 
ferent observables converge at different rates) [ 
[S]. In the NCFC case, one extrapolates to the 
continuum limit as illustrated in Fig. [TJ where 
we show results for the ground state of 12 C as a 
function of N max obtained with a realistic NN in- 
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teraction, JISP16 [E]. The smooth curves por- 
tray fits that achieve the desired independence 
of N max and Q so as to yield the extrapolated 
ground state energy. Our assessed uncertainty 
in the extrapolant is about 2 MeV and there is 
rather good agreement with experiment within 
that uncertainty. The largest cases presented in 
Fig. [T] correspond to N max — 8, where the ma- 
trix reaches a basis dimension near 600 million. 
Computation for the N max = 10 matrix case is 
in progress. Large scale calculations like these 
are performed on leadership-class parallel com- 
puters, at Argonne National Laboratory and at 
Oak Ridge National Laboratory, to solve for the 
low-lying eigenstates and eigenvectors as well as 
to carry out evaluation of a suite of experimental 
observables. For example, one can now obtain the 
low-lying solutions for A = 14 systems with ma- 
trices of dimension one to three billion on 8000 to 
50000 processors within a few hours of wallclock 
time. Since the techniques are evolving rapidly [ 
112] and the computers are growing dramatically, 
much larger matrices are within reach. 



3. Basis Light Front Quantized approach 

To take full advantage of the similarities with 
non-relativistic quantum many-body theory in 
what we will term a "Basis Light Front Quan- 
tized (BLFQ)" approach, we adopt a light-front 
single-particle basis space consisting of the 2-D 
harmonic oscillator for the transverse modes (ra- 
dial coordinate p and polar angle <j>) and a dis- 
cretized momentum space basis for the longitudi- 
nal modes. Adoption of this basis is also consis- 
tent with recent developments in AdS/CFT cor- 
respondence with qcd [ us na na us ■ we 

define our light-front coordinates as x^ — x°±x 3 , 
x 1 - = (a; 1 ,! 2 ) and coordinate pair (p,cf>) are the 
usual cylindrical coordinates in (x , x ). The 
variable x + is light-front time and x~ is the lon- 
gitudinal coordinate. We adopt x + = 0, the "null 
plane" , for our quantization surface. 

The 2-D oscillator states are characterized by 
their principal quantum number n, orbital quan- 
tum number m and harmonic oscillator energy 
f2. It is also convenient to interpret the 2-D os- 
cillator as a function of the dimensionless radial 



length scale variable y/MoHp where Mo has units 
of mass and p is the conventional radial variable 
in units of length. The properly orthonormal- 
ized wavef unctions, & n ,m(p, 4 1 ) = (p<P\nm) = 
fn,m{p)Xm(4>), are given in terms of the Gener- 
alized Laguerre Polynomials, L„ (Mq O p 2 ), by 



f n , m {p) = y/2M Q 



Xm{4>) 




-MoOp 2 /2 



L^(M n P 2 ) 



(i) 



with eigenvalues E n m = (2n + |m| + In 
this 2-D oscillator basis the Fourier transformed 
wavefunctions have the same analytic structure in 
both coordinate and momentum space, a feature 
reminiscent of a plane- wave basis. 

The longitudinal modes, ipj , in our basis are de- 
fined for —L < x ~ < L with both periodic bound- 
ary conditions (PBC) and antiperiodic boundary 
conditions (APBC): 

r/> k {x-) = -Le l £ fc *~, (2) 



where k = 1, 2, 3, ... for PBC (neglecting the zero 
mode) and k = |, §,§>••• in E q-@ for APBC. 
The full 3-D single particle basis state is defined 
by the product form 

^k,n,m{x~, p,4>) = ll)k{x~)$n,m{p,4>)- (3) 

In Figs. [2] and [3] we illustrate some of our basis 
functions. Fig. [2] presents modes for n = 1 of 
the 2-D harmonic oscillator, and Fig. [3] shows the 
transverse mode with n = l,m = joined to- 
gether with the k = ^ longitudinal APBC mode 
of Eq.© and display slices of the real part of 
this 3-D basis function at selected longitudinal 
coordinates, x~ . As one increases the orbital 
quantum number m, pairs of maxima and min- 
ima populate the angular dependence of the ba- 
sis function. Also, as one increases the princi- 
pal quantum number n, additional radial nodes 
appear. As explained in [ HOj and references 
therein, our choice of basis functions is supported 
by the the phenomenological success of the "soft- 
wall" AdS/QCD model which usc s a har- 
monic oscillator potential in the fifth dimension 




of Anti-de Sitter space to simulate color confine- 
ment. The AdS/QCD model, together with light- 
front holography [ [TI3 EE E] > provides a semi- 
classical first approximation to strongly coupled 
QCD. The BLFQ approach in this paper provides 
a natural extension of the AdS/QCD light-front 
wavefunctions to multiquark and multi-gluonic 
Fock states, thus allowing for particle creation 
and absorption. By setting up and diagonalizing 
the light-front QCD Hamiltonian on this basis, we 
incorporate higher order corrections correspond- 
ing to the full QCD theory, and we hope to gain 
insights into the success of the AdS/QCD model. 

4. Cavity mode light-front field theory 
without interactions 

For a first application of the BLFQ approach, 
we consider a non-interacting QED system con- 
fined to a transverse harmonic trap or cavity. The 
basis functions are matched to the trap so the 
length scale of the 2-D harmonic oscillator ba- 
sis is fixed by the trap and finite modes in the 
longitudinal direction are taken with PBC from 
Eq.([2|) (we omit the zero mode). Since we are 
ultimately interested in the self-bound states of 
the system, we anticipate adoption of the NCSM 
method [ 5] for factorizing the eigensolutions into 
simple products of intrinsic and total momentum 



Figure 3. Transverse sections of the real part of 
a 3-D basis function involving a 2-D harmonic 
oscillator and a longitudinal mode of Eq. ([2]) with 
antiperiodic boundary conditions (APBC). The 
quantum numbers for this basis function are given 
in the caption. The basis function is shown for the 
full range — L < x~ < L. Figure from Ref . [ [TO] . 



solutions in the transverse direction. That is, 
with a suitable transverse momentum constraint 
such as a large positive Lagrange multiplier times 
the 2-D harmonic oscillator Hamiltonian acting 
on the total transverse coordinates, the low-lying 
physical solutions will all have the same expecta- 
tion value of total transverse momentum squared. 
Therefore, following Ref.[[T] we introduce the to- 
tal invariant mass-squared M 2 for these low-lying 
physical states in terms of a Hamiltonian H times 
a dimensionless integer for the total light front 
momentum K 

M 2 + P x Px M 2 + const = P + P~ = KH (4) 

where we absorb the constant into M 2 . The non- 
interacting Hamiltonian Hq for this system is de- 
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fined by the sum of the occupied modes i in each 
many-parton state with the scale set by the com- 
bined constant A 2 = 2MgO: 

H = 2M P- 

_ 2M n ^ 2n, + | rTi, | + 1 + 77i 2 /(2A/ r2) 

(5) 

where ttt^ is the mass of the parton i. 

We adopt symmetry constraints and two cut- 
offs for our many-parton states. For symmetries, 
we fix the total charge Z 7 the total azimuthal 
quantum number M t , and the total spin projec- 
tion S along the x~ direction. For cutoffs, we 
select the total light-front momentum, K, and 
the maximum total quanta allowed in the trans- 
verse mode of each many-parton state, N max - 
The chosen symmetries and cutoffs are expressed 
in terms of sums over the quantum numbers of 
the single-parton degrees of freedom contained in 
each many-parton state of the system in the fol- 
lowing way: 

E* = z ( 6 ) 

i 

Y,mi=M t (7) 

i 

E s * = 5 ( 8 ) 

i 

i i 

Ys 2n i + \ m i\ + 1 < N ™ a x (10) 

i 

where, for example, ki is the integer that de- 
fines the PBC longitudinal modes of Eq.Q for 
the i th parton. The range of the number of 
fermion-antifermion pairs and bosons is limited 
by the cutoffs in the modes (K and N max ). Since 
each parton carries at least one unit of longitu- 
dinal momentum, the basis is limited to K par- 
tons. Furthermore, since each parton carries at 
least one oscillator quanta for transverse motion, 
the basis is also limited to N max partons. Thus 
the combined limit on the number of partons is 
vaxa.(K, N max ). Only the case with simultaneous 
increase in both of these cutoffs keeps the problem 



physically interesting at higher excitations since 
this is the only case with unlimited number of 
partons as both cutoffs go to infinity. In princi- 
ple, one may elect to further truncate the many- 
parton basis by limiting the number of fermion- 
antifermion pairs and/or the number of bosons 
but we have not elected to do so here. 

In a fully interacting application, the actual 
choice of symmetry constraints will depend on 
those dictated by the Hamiltonian. For exam- 
ple, with QCD we need to add color and flavor 
attributes to the single particle states and ap- 
ply additional symmetries such as requiring all 
many-parton states to be global color singlets as 
discussed below. Another example is the choice 
to conserve total M t + S rather than conserving 
each separately as chosen here. It is straightfor- 
ward, but sometimes computationally challeng- 
ing, to modify the symmetries in a basis function 
approach such as we adopt here. However, in or- 
der to approach the continuum limit (all cutoffs 
are removed) as closely as possible with limited 
computational resources, one works to implement 
as many of the known symmetries as possible. 

In FigH] we illustrate how the BLFQ basis- 
space dimensions rapidly increase as a function 
of dimensionless state energy for non-interacting 
QED with massless partons. The states are 
grouped to form a histogram according to their 
energy calculated from the chosen Hamiltonian 
in Eq.([5]) where we omit the constant preceding 
the summation for simplicity. We set K = N max , 
and increase them simultaneously. For simplicity 
we consider a case with no net charge Z = 0, 
i.e. for zero lepton number. Thus the cavity 
is populated by many-parton states consisting 
of fermion-antifermion pairs and photons. The 
chosen symmetries are M = and S = 0. 
As seen from Eqs. (|5|9j ). fixing K but increasing 
Nmax leads to a situation where partons carry- 
ing large longitudinal momentum fraction become 
rare, and thus the low-lying modes get eventually 
maximally populated, leading to saturation. As 
shown in Ref. [HQ], the energy at which this sat- 
uration occurs, increases with N max . In the case 
shown here, there is no saturation in state den- 
sity at low energy. With increasing cutoff, there 
is a rapid growth in the number of basis states 
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within each Fock space sector. Overall, there is 
approximately a factor of 20 increase in the to- 
tal many-parton basis states with each increase 
of 2 units in the cutoff. For the non-interacting 
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Figure 4. State density as a function of di- 
mensionless state energy E from BLFQ for non- 
interacting QED in a trap with no net charge and 
for K = N max . The dimensions of the resulting 
matrices are presented in the legend. The states 
are binned in groups of 5 units of energy (quanta) 
where each parton carries energy equal to its 2- 
D oscillator quanta (2n; + |m.;| + 1) divided by 
its light- front momentum fraction (xi = ki/K). 
Figure from Ref . [ HO] . 



Hamiltonian the scale is available through an 
overall factor A 2 = 2M il as described above. 
Without interactions and the associated renor- 
malization program, one cannot relate the scales 
at one set of (K,N max ) values to another. Ulti- 
mately, one expects saturation will arise with in- 
teraction/renormalization physics included as one 
increases the set of (K 7 N max ) values. 

These state densities could serve as input to 
model the statistical mechanics of the system 
treated in the microcanonical ensemble. Of 
course, interactions must be added to make the 



model realistic at low temperatures where corre- 
lations are important. After turning on the inter- 
actions, the challenge will be to evaluate observ- 
ables and demonstrate convergence with respect 
to the cutoffs (N max and K). Independence of 
the basis scale, il, must also be obtained as N max 
increases Ref. [ 18 j. These are the standard chal- 
lenges of taking the continuum limit. 

The approach can be extended to QCD by im- 
plementing the SU(3) color degree of freedom for 
each parton - 3 colors for each fermion and 8 
for each boson. For simplicity, we restrict the 
present discussion to the situation where iden- 
tical fcrmions occupy distinct space-spin single- 
particle modes. We consider two versions of im- 
plementing the global color-singlet constraint for 
the restricted situation under discussion here. In 
both cases we enumerate the color space states 
to integrate with each space-spin state of the cor- 
responding partonic character. In the first case, 
we follow Ref. [[19] by enumerating parton states 
with all possible values of SU(3) color. Thus each 
space-spin fermion state goes over to three space- 
spin-color states. Similarly, each space-spin bo- 
son state generates a multiplicity of eight states 
when SU(3) color is included. We then con- 
struct all many-parton states having zero color 
projection. Within this basis one will have both 
global color singlet and color non-singlet states. 
The global color-singlet states are then isolated 
by adding a Lagrange multiplier term in many- 
parton color space to the Hamiltonian so that 
the unphysical color non-singlet states are pushed 
higher in the spectrum away from the physical 
color single states. To evaluate the increase in 
basis space dimension arising from this treatment 
of color, we enumerate the resulting color-singlet 
projected color space states and display the re- 
sults as the upper curves in Fig. [5] 

In the second case, we restrict the basis space to 
global color singlets and this results in the lower 
curves in Fig. [5] The second method produces 
a typical factor of 30-40 lower multiplicity at the 
upper ends of these curves at the cost of increased 
computation time for matrix elements of the in- 
teracting Hamiltonian. That is, each interacting 
matrix element in the global color-singlet basis is 
a transformation of a submatrix in the zero color 
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projection basis. Either implementation dramat- 
ically increases the state density over the case of 
QED, but the use of a global color-singlet con- 
straint is clearly more effective in minimizing the 
explosion in basis space states. 
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Figure 5. Number of color space states that ap- 
ply to each space-spin configuration of selected 
multi-parton states for two methods of enumer- 
ating the color basis states. The upper curves are 
counts of all color configurations with zero color 
projection. The lower curves are counts of global 
color singlets. Figures from Ref. [ [10] . 



5. Elements of the interacting theory 

In the interacting theory a primary concern will 
be to manage the divergent structure of the the- 
ory. There are two possible locations for diver- 
gences in a Hamiltonian basis function approach: 



(1) the matrix elements themselves diverge, or 

(2) the eigenvalues diverge as one or more cut- 
offs are removed. In our planned cavity field the- 
ory application with interactions, we will manage 
these divergences with the help of suitable coun- 
terterms and boundary conditions. The infrared 
divergences in light-front momentum arising in 
both the fermion to fermion-boson vertex as well 
as in the instantaneous fermion-boson interaction 
are expected to be well-managed by previously 
defined counterterms [ [3] suitably transcribed for 
the transverse basis functions we have adopted. 
We anticipate this prescription will work since 
the longitudinal modes we adopt are similar to 
those used in Ref. [ [3] . Finally, we expect the 
transverse ultraviolet divergences to be suitably- 
managed with our basis function selection. As an 
alternative scheme for comparison we also plan 
to adopt a second approach that involves a re- 
cently proposed sector-dependent coupling con- 
stant renormalization scheme [ [20] . Another al- 
ternative which we may adopt uses the Pauli- 
Villars regulator [|2T]. 

Since we are introducing a basis-function ap- 
proach for the transverse degrees of freedom, we 
need to investigate convergence rates with in- 
creased cutoff of the transverse modes N max . 
Here, in a simple set of examples, we outline how 
we can search for additional sources of divergence 
with the help of a perturbation theory analysis. 
For a first investigation, we have examined the be- 
havior of various sets of matrix elements for the 
fermion to fermion-boson vertex. For the pur- 
pose of this investigation, we adopt periodic (an- 
tipcriodic) boundary conditions for the longitudi- 
nal modes of the bosons (fermions) and we hold 
spin projections fixed for initial and final states. 
We then adopt specific values for the longitudi- 
nal momentum fractions observing the conserva- 
tion rule. The trends we examine should not be 
sensitive to the specific values adopted. 

Consider the second order energy shift, AE, in- 
duced on a single parton in the transverse mode 
(n, m) by its coupling V to partons in higher en- 
ergy transverse modes (n',m'): 

f \(n,m\V\n',m')\ 2 , , , . , 
A£„, m « / IV ' 1 ' ' ^p{n')dn\ (11) 
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where 

n' = 2n' + \m'\ (12) 

Eni, m , = (n' + l)(l (13) 

p(n') = n' + l (14) 

and the degeneracy is taken into account in p. 

Thus, according to perturbation theory, we ex- 
pect a UV divergence if the matrix element falls 
off too slowly with increasing n! . In particular, if 
the matrix element falls approximately as (n')~^ ; 
then we expect a logarithmic divergence since the 
integrand will have a net (n')~ 1 dependence. If 
the falloff is even slower then we encounter a more 
serious divergence. Another possible source of a 
log divergence could arise within the selected sum 
over ml in which case p, the level density factor 
in the integrand, is unity. Then, if the matrix el- 
ements for fixed n, n' are approximately constant 
with increasing m', we again find a log divergence 
in the sum over m! . 

We portray in FigsE] and [7] representative se- 
quences of how these off-diagonal matrix ele- 
ments behave as one increases the difference in 
the initial and final state principal quantum num- 
bers. We also portray two interesting cases where 
the fcrmion and fermion-boson principal quantum 
numbers track each other, cases that do not en- 
ter a perturbative analysis. Note that we have 
limited the illustrations to the transverse com- 
ponents of our matrix elements. We also select 
cases where the fermion spin is flipped, cases that 
are proportional to the fermion mass. We set 
the fermion mass to unity so the results are ex- 
pressed in units of the fermion mass. We fur- 
ther limit our presentation to the case where all 
partons remain in the orbital projection quan- 
tum number zero state and the 2-parton (fermion- 
boson) states have each parton in the same trans- 
verse state. For the non-perturbative illustrative 
cases, we display the matrix element trends where 
all partons remain in the same transverse mode, 
(n,0). 

When the single fermion state has an even value 
of the principal quantum number n as shown in 
Fig. [SI the matrix elements appear to be well- 
behaved either when the 2-parton configuration 
is the lowest accessible case ((0000|) or when the 
each of the two partons resides in the state with 



the same principal quantum number as the single 
fermion state. We demonstrate anticipated good 
convergence with increasing n by showing that 
the matrix elements, when multiplied by y/n + 1, 
still fall with increasing n. For the case when the 
single fermion state has an odd value of the prin- 
cipal quantum number n as shown in Fig. [7j the 
situation is somewhat different. For the matrix 
element set entering a perturbative analysis, the 
matrix elements fall to zero with increasing n suf- 
ficiently fast that multiplying by y/n+1 does not 
significantly distort the trend to zero. However, 
the large n behavior of the fermion-boson matrix 
element, with all partons at the same n, is seen 
to go approximately as y/n + 1. This is best seen 
in Fig. [7] where the matrix elements are mul- 
tiplied by y/n+1 and the result appears to be a 
nonzero constant at large n. Since this trend does 
not appear in a second order perturbation theory 
analysis, we must await the full Hamitonian di- 
agonalization to better understand its role in the 
convergence with increasing N max . To date, 
we have found no matrix element trends with in- 
creasing transverse energy that would imply new 
divergences. All such matrix elements sets we ex- 
amined, within a perturbative analysis, fall faster 
than 1 / y/nJ . Also, we found no sets that remained 
constant with increasing m' (and thus m) holding 
n, n' fixed. As a result of this initial analysis, we 
anticipate that straightforward adoption of coun- 
terterm methods previously introduced [[3] will be 
sufficient for managing the identified divergences 
in BLFQ. 

As another example of the interacting cavity 
mode QED, we present in FiglS] the eigenvalue 
spectra for a Hamiltonian matrix that includes 
lepton and lepton-photon Fock spaces. The lep- 
tons are now massive with fhi — 0.511 MeV, 
and the HO parameters are f2 = 0.1 MeV and 
Mq = fhi . The other basis parameters are chosen 
as K = 3, N max = 2 and Mj = 1/2. The eigen- 
values are shown as a function of the strength 
of the coupling with respect to the coupling at 
QED, with <7qed = y/4^ra. Since the instanta- 
neous fermion interaction term only contributes 
when all the leptons and photons of the process 
have their spins aligned, this particular case only 
involves the photon emission vertex interaction 
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Figure 6. Behavior of representative fermion to 
fermion-boson matrix elements in BLFQ. The 
quantum numbers specifying the parton trans- 
verse modes (rii,mi) in the matrix elements are 
given in the legend. Only the transverse mode 
contributions to the matrix elements are shown. 
Results are also shown with a multiplicative fac- 
tor of y/n + 1 applied to help search for a loga- 
rithmic divergence by obtaining a resulting flat 
behavior, when it occurs. Overall matrix element 
normalization depends on the specific values of 
light-front momentum fractions carried by the in- 
teracting partons. Figure from Ref . [ [TU] . 



(see e.g. Ref.[[22]). Application of the sector- 
dependent renormalization [ [20] will shift these 
eigenvalues such that the lowest one will obtain 
the value E — m ^ 8 " , where to is the mass 
of the lowest state fermion. As we increase the 
cutoffs K,N max , the evaluation of the electron 
anomalous magnetic moment will serve as an im- 
portant precision test for our method [ [22] . 

6. Conclusions 

Following successful methods of ab initio nu- 
clear many-body theory, we have introduced a ba- 
sis light-front quantization (BLFQ) approach to 
Hamiltonian quantum field theory. With a cavity 
mode treatment of massless non-interacting QED 



Figure 7. As in Fig[6] Note that for one of 
the cases shown here, the resulting matrix ele- 
ments vanish with increasing n while the other 
case shows a constant trend but does not enter a 
second-order perturbative analysis. Figure from 
Ref. [HO]. 



we showed how the number of many-parton ba- 
sis states exhibit a dramatic rise as the cutoffs 
are elevated. In order to extend our method to 
QCD, we presented two methods for treating the 
color degree of freedom and demonstrated sample 
measures of the efficiency gains for global color 
singlets over basis states with color-singlet pro- 
jection alone. For the interacting cavity mode 
QED we outlined our approach to managing the 
expected divergences in a manner that will pre- 
serve all the symmetries of the theory. An ini- 
tial inspection of the interaction vertices of QED 
in the BLFQ approach shows smooth behaviors 
that, following a second-order perturbative anal- 
ysis, are not expected to lead to divergences. It 
appears that the cavity mode treatment, with 
the type of basis spaces we have selected, will 
encounter the divergences in a more subtle fash- 
ion as cutoffs are elevated. The evaluation of the 
electron anomalous magnetic moment will be an 
important test to our method as we expand our 
basis space. 
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Figure 8. Eigenvalues of a cavity mode QED 
Hamiltonian, which includes lepton and lepton- 
photon Fock spaces, as a function of the strength 
of the coupling (<7qed = %/47ra). The HO param- 
eters are ft = 0.1 MeV and M = 0.511 MeV, and 
the other basis parameters are K — 3, N max = 2 
and Mj = 1/2. 
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